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Abstract 

We justify the practical use of the Shuvaev integral transform approach to calculate 
the skewed distributions, needed to describe diffractive processes, directly from the con- 
ventional diagonal global parton distributions. We address doubts which have been raised 
about this procedure. We emphasise that the approach, on the one hand, satisfies all the- 
oretical reqirements, and, on the other hand, is consistent with DVCS data at NLO. We 
construct an easily accessible package for the computation of these skewed distributions. 



1 Motivation 

Skewed parton distributions are needed to calculate exclusive diffractive production, such as 
7 + p — > VM + p for light or heavy vector mesons, see e.g. [H [21 El S], central exclusive 
diffractive Higgs boson production at the LHC [5], etc. For all these diffractive processes we 
need skewed parton distributions at small values of £ <C 1, and small to medium scales. Data 
for diffractive J/ip production, for example, test the gluon in the range £ ~ x = 10~ 4 . . . 5 • 10 -3 




Figure 1: A schematic diagram showing the variables for the off-diagonal parton distribution 
H(x,£) where x\^ — x ± £. 



and effective scales fi 2 ~ 2... 8 GeV 2 [6]. For exclusive Higgs production at the LHC, the 
regime of £ ~ x ~ 10~ 2 and scales of the order 5 GeV 2 will be relevant, see [5]. 

Moreover it has been proved, using dispersion relations [7], that only distributions in the 
space-like region \x\ > £ are needed to describe these processes. The variables are defined in 

Fig.m 

There are insufficient experimental data to determine the skewed parton distributions with 
an accuracy comparable to that of the global parton analysis of the conventional (diagonal) 
distributions. Fortunately, for small £, the skewed distributions can be computed accurately 
just from the knowledge of the known integrated conventional distributions. At first sight such 
a simplification looks surprising. On the other hand, we know that the anomalous dimensions 
which describe the evolution of the Gegenbauer moments, Gn, of the skewed distributions 
are equal to the corresponding anomalous dimensions of the conventional Mellin moments, 
M N [H [9] . This is a consequence of conformal invariance of the evolution equations. Strictly 
speaking, conformal invariance is only valid at leading order (LO). Already at NLO [TUl [TTj . it 
is violated by the running of as, leading to a mixture of the operators at adjacent orders. But 
let us start with LO. 

At LO we have equality of the anomalous dimensions of Gn and Mn- Moreover, due to the 
polynomial property [121 [131 [H]i 

N 

Gn = £ c ^ 2n > C 1 ) 

n=0 

we have c$ = M^. That is, from the conventional global parton analyses we can determine 
all the Gegenbauer moments of the skewed parton distributions at small £ with an accuracy 
of 0(£ 2 ). Then it is simply a technical problem to calculate the x distribution of the skewed 
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partons from the known moments. This mathematical problem is solved by the Shuvaev trans- 
form [T5] (which we give explicitly later). Note that at £ = 0, in the diagonal case, there is 
no mixture of the Mellin moments during the evolution. Thus for the skewed distributions 
such mixing must also vanish at £ = 0. Indeed the mixture of the different operators, due to 
violation of conformal invariancqj, is of 0(a s £) [HE \TT\ [T6| ITT] . Thus the Shuvaev transform 
can be used at NLO with accuracy 0(£), which is sufficient for the description of all diffractive 
processes of interest. 

This procedure has been called into doubt [18]. First, we discuss the reason for the doubt 
and then explain why it is not a problem in practice. The apparent problem is that to obtain 
the x distributions from the moments Gn we must analytically continue the moments into the 
complex iV-plane. If there is a singularity in the right-half plane, then it may generate a non- 
negligible correction of 0(£/x), instead of the 0(£ 2 ) correction which came from the difference 
between the Gegenbauer and Mellin moments. This large 0(£/x) correction would destroy the 
practical use of the Shuvaev transform. So to justify the Shuvaev transform we require the 
absence of singularities in the right-half N plane. 

Now, an arbitrary singularity in the right-half plane will violate the polynomial property, 
and so at first sight the danger is removed. However, it has been shown by Radyushkin [T9l |20~1 
[2TJ, [22] that it is possible to form 'double distributions' which are not identical to the Shuvaev 
transform, but which still satisfy the polynomial property. The non-polynomial contributions, 
generated by different singularities of these double distributions in the right-half plane, com- 
pensate each other to guarantee polynomiality. So the danger remains. On the other hand, 
there is no singularity in the right-half plane in the anomalous dimensions which describe the 
q 2 evolution of the Gegenbauer moments. Hence the extra singularities must come only from 
the input distribution. Now, it is natural to describe the input distribution at low x in terms 
of the Regge approach, which is successful in the description of high energy interactions at low 
scales where the conventional (collinear) DGLAP evolution for the skewed distributions starts. 
In the Regge approach there are no singularities in the right-half plane (J > 1) in the space-like 
{x > £) domain. 

Indeed, in the Regge limit we are concerned with the leading log(l/x) summation. Hence, 
in Fig. [2] we have the strong ordering 

%i 3> x' i+1 , that is \xi — x[\ = 2£ <C Xi for i < n. (2) 

The inequality £ <C X, means that the lower part of the diagram is described by the diagonal 
distribution. The only £ dependence comes from the uppermost x n cell. All other cells have 
Xi ^> £. Now, the uppermost cell already satisfies conformal invariance. Moreover, at NLO level, 
this uppermost cell is calculated exactly in terms of the NLO coefficient functions. Therefore 
there is no opportunity for a new singularity in the right-half j > 1 plane. 

This violation occurs due to the dependence of a s on the dimcnsionful parameter Aqcd- It is essentially 
trivial. In principle, it could be accounted for by using the 'correct' argument of the QCD coupling in the 
uppermost cell of Fig. [TJ namely a s (k^) and not a s (Q 2 ). 
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Figure 2: Schematic diagram describing deeply virtual Compton scattering (DVCS) in terms of 
generalised parton distributions. The skewing (or £ dependence) originates from the uppermost 
cell; the part of the diagram enclosed by the dashed line has diagonal form due to the strong 
ordering of x in the Regge limit. This diagram makes clear why, in the low x region, gener- 
alised parton distributions may be computed directly from the well known 'global' diagonal 
distributions. 



With the above physically motivated conjecture, the Shuvaev transform survives to be of 
practical value. We emphasise that the conjecture, of the absence of additional singularities in 
the right-half plane, does not follow from first principles or symmetry arguments, but rather 
from the assumption of the Regge form of the low x input. 

Note that in the time-like region, that is at \x\ < £, we have such singularities. They 
correspond to the wave functions of resonances with large spin. Indeed one may add to the 
distribution, given by the Shuvaev transform, an additional function with support only on the 
interval \x\ < £. However there is no such contribution in the calculation of the skewed partons 
via the Shuvaev transform in the space-like domainal 

Another argument in favour of the Shuvaev prescription is the good fit of the DVCS HERA 
data obtained at NLO and NNLO in [23J. This fit was based on the approximation that there 
is only one pole. Recall also the model [21] based on the Shuvaev transform and the successful 
description of the DVCS HERA data in the dipole model [25J where the ratio, R g , of the skewed 
to diagonal gluon distribution given by the Shuvaev transform was used. 

2 As was pointed out in [24j . such 'meson- wave- function- like' contributions lead to extra terms in the inverse 
Shuvaev transform, which we do not use here. Also note that similar extra terms in the effective diagonal 
functions as proposed in |24j are irrelevant for our case. 
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In [26] the Shuvaev transform was criticised to be "unpractical" or "too complicated" , but 
see also the original discussion in [27|. Here we present and describe a package that allows 
for the simple computation of the skewed distributions at small x, £ based on the Shuvaev 
transform. 



2 Description of the Shuvaev transform 



The generalised parton distributions (GPDs) are denoted [12][T3l[14] by Hi(x,£), i = q,g where 
£ is the skewing parameter, with x 1>2 = x±£ as defined in Fig.[TJ In fact H(x,£) = H(x, £, t, fi 2 ) 
for partons emitted and absorbed at scale fi 2 , with the momentum transfer t = (p — p') 2 , and 
— 1 < x < 1. The values of t and £ do not change as we evolve in /i 2 . In the limit £ — > the 
skewed distributions will reduce to the usual diagonal partons: 



H q (x,0) 
H g (x,0) 



q{x) 



x g(x). 



for x > 
-x) for x < 



(3) 



We briefly discuss the double-integral (that is, the Shuvaev transform) used to calculate the 
skewed distributions for £ <C 1 in the space-like region \x\ > £. First we introduce an auxiliary 
function f^(x,t), whose Mellin moments are equal up to normalisation to the Gegenbauer 
moments of the skewed function. The function formed from the DGLAP evolution of the 
auxiliary function is referred to as the 'effective diagonal function', f(x'). Neglecting the t 
dependence, we must find the kernel K(x, £; x') relating f(x') to the skewed distribution, 



H(x,0 = J dx'K(x,Z;x')f(x'). 



(4) 



It is convenient to first weaken the singularity in the x' integral by an integration by parts. 
With this it can be shown that [28] 



L dx ' 



Im 



7T 



ds 



o y(s) y/l-y{s)x' 



d / q(x') 



dx' 



2 ^ r 1 ds(x + gi-2s)) 
vr Jo y(s) y/l-y(s)x' 



dx' \ \x'\ 
d fg{x') 



dx' V \x 



4s(l-s) 



(5) 



(6) 



(7) 



Equations © and © are (up to a variable substitution z = l/(x'y(s))) the form used to 
compute the Shuvaev transform in the package decribed below. 



5 



Incidentally, it is also possible to solve the inverse problem. That is to obtain the diagonal 
distribution from a known skewed distribution H(x,£) at a given value of £. Although this is 
not needed for our discussion, for completeness, we give the details in the Appendix. 

Note that the transforms, fl5]) and flSD, say nothing about the t (or p T ) dependence of the 
GPDs. Strictly speaking they are written for pt = 0. Usually the factorisation 

H(x^; PT ) = H{x,0F{pr) (8) 

is assumed, where F(pt) is just the proton form factor. 



3 Approximate predictions of GPDs for small x and £ 

We see that (jSJ) and (jSJ) determine the behaviour of the skewed distributions in the small x, £ 
domain entirely in terms of the diagonal distributions. Before we perform the exact evaluation of 
these integral expressions for the GPDs, it is informative to recall that approximate expressions 
can be obtained by making the physically reasonable small x assumption that the diagonal 
partons are given by 

xq(x) = N g x~~ Xq , xg(x) = N 9 x~ Xb . (9) 
Then we can perform the x' integration analytically^ We obtf 



4s(l 



tarn 

) l Ai+1 
x + £{l-2s)_ 

with i = q or g, and where p = and 1 for quarks and the gluon respectively. 



G(t) (10) 



At first sight it appears that for singlet quarks (where X q > and p = 0) we face a 
strong singularity in integral ffTUl) when the term D = x + £(1 — 2s) —>■ in the denominator. 
Fortunately the singlet quark distribution is antisymmetric in x. To obtain the imaginary part 
of the integral ([SD we must choose x' > for D > and x' < for D < 0. Therefore we 
must treat ffTU]) as a principal value integral and take the difference between the D —>■ 0+ and 
D — ► 0— limits. Thus the main singularity is cancelled and (|T0|) becomes integrable for any 
A, < 1. 

Note that the dominant contribution to the x' integrations of (jSJ) and ([6]) comes from the 
region of small x' ~ x, £. Indeed with the input given by ([9~lh the integral for the quark 



3 With the substitution z = X/x'y(s) one can then use 



1 , a+3 ,„ , x r(A + |) r(£ 







r(A + 3) 



where we have set the lower limit to zero. We have checked numerically that this is a very good approximation, 
with accuracy of the order of 10~ 4 for A ~ 0.2 in the small x region. (That this approximation is good for small 
£ is in line with the findings of [24].) 
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Figure 3: The off-diagonal to diagonal ratio, R a , at small x = £ versus the power A which 
specifies the x~ x behaviour of the input diagonal parton as in ([9]). Note that the quark singlet 
ratio has been divided by two. 



distribution has a strong singularity at small x' 

i q ~ / dxvr A - 3 im i ; i . (ii) 

However when we take the imaginary part, the x' integration is cut-off by the theta function 
8(x' - l/y{s)) at 

x' = l/y(s) ~ x + Z(l-2s), (12) 

which implies that x' is always greater than x/2. So we obtain the small £ behaviour I q ~ ^ _Ag_1 , 
and the distribution has the form 

H q (x,£) = C**- 1 F q (x/t). (13) 

Similarly it follows that H g = £~ Xs F g (x/£). 

It is illuminating to evaluate the ratio R of the GPD to its diagonal parent distribution; 
that is 

H(x,0 

R - or (14) 

where the only free parameter is A, the exponent which fixes the x~ x behaviour of the input 
diagonal partons, as in Notice that on account of f[T3"j) the ratios R at small x and £ are a 
function of only the ratio of the variables 

Assuming the pure power behaviour of the partons, the ratios at x = £ are given explicitly 
in analytic form as 

Ra m0 2^ T(A + 5/2) 

H(2£,0) v^F r(A + 3 + p)' 1 ] 

where p = for quarks and p — 1 for gluons. These ratios are plotted in Fig. [3]as a function of A. 
We see that the off-diagonal or 'skewed' effect (the ratio R a ) is much stronger for singlet quarks 



than for gluons. The explanation is straightforward. At low x the distributions are driven by 
the double leading logarithmic evolution of the gluon distribution. At each step of the evolution 
the momentum fractions x^ are strongly ordered (x[ ^> Xx,x' 2 3> x 2 on Fig. 1). For gluons it 
is just the "last splitting function" P gg (x2, x' 2 ; £) which generates the main £ dependence, or 
skewedness, of the distribution. However for the sea or singlet quarks it is necessary to produce 
a quark with the help of P qg at the last splitting. The splitting function P qg has no logarithmic 
1/z — x 2 /x 2 singularity and so x 2 is the order of x 2 . Consequently both the splitting functions 
P qg (x2, x' 2 ',C) and P gg (x' 2 , x 2 ; £) generate the asymmetry of the off-diagonal distribution. Hence, 
at low x, the singlet quark has a much stronger off-diagonal effect than the gluon. 

We emphasise that the analytical expression for the ratio R a , ( fl5l) . is valid in the limit 
x = £ only and that it assumes a pure power behaviour of the diagonal partons, (JHJ). As will be 
discussed below, the latter is a good approximation for global fit partons in the small x regime. 
However, for x > £ the difference between the result of the complete (double integral) Shuvaev 
transform, (J5], [6]) , and the approximate analytic formula for R a , (fT5l) . based on the limit x = £, 
is quite large. This is evident from Fig. [H] below. 



4 Evaluation of the GPDs using the Shuvaev transform 

In the above section we have obtained an approximate determination of the GPDs valid for 
x ~ £. We now perform a precise evaluation of the Shuvaev transforms of ([5]) and ([6]), valid 
for all small x, £, and compare our results to the approximation based on (jTUl) . This is done 
in Fig. H] for the gluon and in Fig. [5] for the sums of u, d, s quarks and antiquarks, using 
MSTW2008NLO (conventional global fit) parton distributions as input. Thus our GPDs 
correspond to using the MS renormalisation scheme and the General Mass Variable Flavour 
Number Scheme adopted by MSTW. 

In the lower panels of the figures the effective powers A are plotted as a function of x for 
three different scales, /i 2 = 2.5, 10, 50 GeV 2 . In the upper panels both the ratios 

H(x,0) K ' 

for the full Shuvaev transform and the corresponding analytical approximations 
~ H(x/2,x/2) 2 2A + 3 T(A + 5/2) . . , x 

R a = \J ' ; = — \ — '— with A evaluated at x) 17 

H(x,0) v^F T(A + 3+p) v ; K J 

are given, again for the three scales. It is clear that, depending on the values of x and the 
scale, the deviation from a pure power can lead to a sizeable difference between the analytic 
approximation and the full result. However, for small x < 2 ■ 10 -3 and not too small scales 
the deviation is quite small in the case of the global fit partons MSTW2008NLO. Actually the 
difference is smaller than it appears at first sight in Fig. [5] since the vertical scale does not 
extend to zero. 
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Figure 4: Lower panel: Values for the effective power A (evaluated at x) of the MSTW2008NLO 
gluon for the three scales /i 2 = 2.5, 10, 50 GeV 2 . Upper panel: Analytic approximation 
(dotted lines) as defined in ffTT|) compared to the ratio R g (solid lines) as defined in ffl6|) for the 
full Shuvaev transform. 
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Figure 5: Lower panel: Values for the effective power A (evaluated at x) for the sum of the 
MSTW2008NLO u, d, s quarks (dashed) and antiquarks (dotted lines) for the three scales /i 2 = 
2.5, 10, 50 GeV 2 . Upper two panels: as defined in (TlTI) shown as dotted lines for quarks (top 
panel) and antiquarks (middle panel) compared to the ratio R q for the full Shuvaev transform 
as defined from (Tl6l) (solid lines). 
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Figure 6: Left panel: Analytic R a g from ( fl5l) (fat dots) compared to the ratio R g from ( JT41) 
(solid lines) using the full Shuvaev transform flS]), at the three scales /i 2 = 2.5, 10, 50 GeV 2 for 
the MSTW2008NLO gluon and £ = 10~ 3 as a function of Right panel: same as left panel 
but for the sum of u, d, s quarks and antiquarks. The insert shows a blow-up of the regime 
x ~ f. 
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Figure 7: Comparison of skewed (solid) and diagonal (dotted lines) gluon (left) and singlet 
quark (right panel) MSTW2008NLO distributions for £ = 10~ 3 and the three scales /i 2 = 
2.5, 10, 50 GeV 2 as a function of x/£. Also shown is the product of the diagonal distributions 
with the analytical skewing enhancement factor R a obtained in the x = £ limit from (1151) 
(dashed lines). 
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Figure 8: Skewed (solid and dash-dotted) and diagonal (dotted lines) MSTW2008NLO (left) 
and CTEQ6.6M (right panel) gluon distributions at scales /z 2 = 2.5, 10, 50GeV 2 . The dashed 
lines show the skewed gluon if the analytic approximation R a from (fTTj) is used (MSTW only). 



Frequently, skewing corrections are taken into account via simple multiplication with the 
ratio R a . One may fear that this approximation of 'maximal skewing' overestimates the real 
skewing effect which comes from integrating over the whole possible range £ < \x\. Possible 
effects are exemplified in Fig. [6] for MSTW2008NLO gluons (left) and the sum of u, d, s quarks 
and antiquarks (right panel). Here the limit R a from (fT5l) (fat dots) is compared to the full 
ratio R from (THj) (solid lines) as a function of x/£ for £ = 10~ 3 at three different scales, 
/z 2 = 2.5,10,50 GeV 2 . 

In Fig. [7] a similar comparison is shown for the skewed gluon (left) and singlet quark distri- 
butions (right panel) as a function of x/£ for £ = 10~ 3 at the three scales. The solid lines show 
the results for the full Shuvaev transforms, H g (x,£) and (x + £)H*(x, £), whereas the dashed 
lines are obtained as the product of the analytical skewing factors R a (in the limit x = £) with 
the diagonal partons evaluated at x + £. The dotted lines show the diagonal MSTW2008NLO 
partons for comparison. It is clear from Figs. [6] and [7] that, at least in the cases under con- 
sideration, the overestimate of skewing effects through simple multiplication with R a could be 
sizeable, although as expected the agreement when x = £ is good. 

In Figs. [8] and [9] we finally show the skewed gluon, H g (x/2,x/2), and singlet quark distri- 
butions, xHg(x/2,x/2) (solid lines), compared to the corresponding diagonal partons (dotted 
lines) as a function of x and at three scales for both MSTW2008NLO [29] and CTEQ6.6M [30]. 
For MSTW, the analytical approximation for the skewing using the factor R a is also shown 
(dashed lines). 

Some features of these results are especially noteworthy. First, we emphasise that for x ^> £ 
the skewed distribution H(x,^) becomes close to the diagonal distribution. This is trivial, since 
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Figure 9: Skewed (solid) and diagonal (dotted lines) quark singlet distributions at scales /z 2 = 
2.5, 10, 50 GeV 2 , for MSTW2008NLO and CTEQ6.6M as indicated in the legend. Shown as 
dashed lines (MSTW only) is the approximation using the analytic skewing factor R a from 
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for x ^> £ we can neglect the £ dependence. Secondly, we see a much larger skewing enhancement 
for the quarks than for the gluons. Indeed, it is known that, in the leading In 1/x approximation 
(LLA), the skewed distributions are equal to the diagonal ones [21] • As mentioned above, this 
arises because, in the LLA, the longitudinal fractions, x, of the momentum are strongly ordered. 
Already at the first evolution step we have x 3> £. Hence the £ dependence becomes negligible. 
Now for gluons, with spin 1, the LLA correponds to a flat xg(x) = constant distribution. With 
a parameterisation of the form xg ~ x~ A , we see from Fig. H]that for A > (A < 0) the skewed 
distribution is enhanced (suppressed) in comparison with the diagonal gluon. On the other 
hand, for a pair of t-channel quarks with spin |, the LLA corresponds to a flat q(x) = constant 
behaviour, see ([3]). That is, for the form xq ~ x~ x , we have no skewed effect if A = — 1, but a 
large enhancement for small A close to zero. In the latter case, with A > 0, the structure of the 
loop integration is such that it prefers to transfer the major part of the momentum flow along 
one quark propagator. 

This behaviour is demonstrated in Figs.[6]and[7], at the lowest scale. We see that the Shuvaev 
transform actually suppresses the gluon. This is expected when we look at the x behaviour 
of the diagonal gluon; inspection of Fig. H] shows that X g becomes negative here. This is also 
demonstrated in the fi 2 = 2.5 GeV 2 curves in Fig. [81 For MSTW gluons, we have a suppression 
of the diagonal gluon, whereas for CTEQ gluons, there is no net effect at this particular choice 
of parameters. This behaviour is demonstrated with the lowest curve for x < 3 • 10~ 2 in the 
lower panel of Fig. [HJ in this regime CTEQ partons have \ g ~ 0. As the scale increases, the 
transform enhances the diagonal gluons. 

We emphasise that the use of the xg, xq oc x~ x forms does not mean that we assume 
power-like asymptotics for x — > 0. Clearly, this oversimplified parameterisation can be valid 
only in a limited range of energy or x. We have used it here simply to qualitatively illustrate 
some of the main features of the skewed effect. Actually, in practice, and in the supplied grids, 
we determine the GPDs directly in terms of the well known diagonal distributions, which have 
much more complicated x structure, using the full Shuvaev transforms, (jSJ) and ([6]). 

4.1 Leading order skewed distributions 

For completeness, we also calculate LO generalised parton distributions. The LO integrated 
distributions have a steeper 1/x behaviour than those at NLO; that is, the values of A are larger. 
Therefore the skewed effect is larger. All other features of the LO generalised distributions are 
qualitatively the same as those at NLO, see Fig. [10] in comparison to Figs. [8] and [9j 

The large numerical difference between the LO and NLO integrated gluon is due, first, to 
the absence of the LO coefficient function for 7*g splitting, C^J = 0; and, second, to a singular 
1/z term in the quark-quark splitting function P qq which is present at NLO, but absent at LO. 
To compensate for these absences at LO, we need to 'artificially' enhance the input LO gluon 
at low x. It is, therefore, not necessarily true that the LO contribution to another process can 
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Figure 10: Skewed (solid and dashed) and diagonal (dotted lines) LO partons from MSTW2008 
as a function of x at scales /z 2 = 2.5, 10, 50GeV 2 . The dashed lines show the skewed partons 
if the analytic approximation R a from (fTTj) is used. Left: gluons, right panel: u, d, s quark 
singlet. 

be effectively replaced by the same enhanced gluons. For this reason we prefer to use NLO 
partons. 

Recall that the description of DVCS data in terms of the LO skewed distribution, given by 
the Shuvaev approach, was not good [32J, while at NLO we observe agreement with the data 
[23J. The data have also been described in terms of a dipole model [25J. Here, the interaction 
of the gluon with the quark-antiquark pair produced by the photon plays the role of the NLO 
coefficient function C ig , while the skewed effect of the NLO gluon was calculated using the 
Shuvaev prescription. Therefore, actually, the dipole approach is close to a NLO treatment of 
DVCS. 

4.2 Grid interpolation package 

As the numerical computation of the full Shuvaev transform requires some care and can be too 
slow for applications, we provide a simple and fast interpolation routine in Fortran77. This 
includes reasonably small (approx. 3.4 MBytes) grid files for the LO and NLO input global 
analyses MSTW2008 [29], MRST2004 [33] and CTEQ6.6 [30]. These contain the information for 
quarks, antiquarks and gluons. The interpolation routine and grid files can be downloaded from 
http://www.maths.liv.ac.uk/TheorPhys/RESEARCH/pubcodesThtmrj0 The grid files use 85 
points in x where 4 ■ 10 -5 < x < 1, together with 43 points in the ratio £/x where < £/x < 1 

4 Grid files for other partons can be created upon request. 
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and 20 points in q 2 where 1.25 GeV 2 < q 2 < 80 GeV 2 . Using £/x as opposed to £ as a parameter 
increases the interpolation accuracy as the region £ < x, where the skewed distributions are 
rising sharply, can be sampled more densely than regions ^ <i, where they are flatter. The x 
grid points were chosen on a log scale, and the transform uses a linear interpolation on log scales 
for the three parameters (apart from the bin which includes £/x = 0, where the interpolation 
is performed on a linear scale). The evolution in q 2 is very smooth, so fewer points in q 2 are 
needed. 

The accuracy of the results from the interpolation on the grid is always better than 0.2% in 
the small x regime (and still better than 1% for x > 10~ 2 ), and typically much better for £ <C x 
or in the case of gluons. Based on convergence tolerances of the double integrals evaluated in 
(J5J [6]), and the numerical derivatives (of the interpolated input partons) of one per mille, we 
estimate the accuracy of the grid points to be better than 0.2% in the small x regime. 

Note that while in the figures above we have shown the skewed partons only up to x = 0.1, 
the grid files contain information up to x = 1. Of course at large x and £ we cannot justify 
the results, while at small £ and large x the skewed distributions are approaching the diagonal 
ones. 

5 Comparison with an alternative approach 

In general there are two possibilities to parameterise GPDs. One is to relate skewed distributions 
to the well known diagonal (global) parton distributions. This can be done with the help of the 
double distributions proposed by Radyushkin [2U1ET]- The general form contains an arbitrary 
new function (up to the normalisation condition), but by construction it reduces to the diagonal 
distributions in the limit £ — > 0. The Shuvaev transform is a particular case of this approach, 
with the advantage that after a physically reasonable assumption (that is, no singularities in the 
right-half j plane) it gives unique GPDs in the low x domain in terms of global diagonal partons, 
without any new parameters. The physical reasonableness of this assumption is evident from 
Fig. [2] and the accompanying discussion. 

An alternative approach is to fit the available data, corresponding to both skewed and 
diagonal distributions, using an ansatz or model for the input GPDs at some starting scale 
Qo- This approach was used in a recent paper [34J, where the moments of the GPDs were 
parameterised in terms of beta (B) functions motivated by an SO (3) partial wave model, see 
Section 3.2 of [31]. In order not to lose the statistical significance of the DVCS data, which 
corresponds to the skewed distributions, only a small subset of global (F 2 ) data were included 
in the analysis. Moreover, we stress that the DVCS data do not cover the whole kinematic 
domain of x and £, but rather correspond to the special situation \x\ = £. It was emphasised in 
[34J that it is impossible to describe HERA data with a model with one leading SO (3) partial 
wave, especially at LO. Based on this observation, the authors claim that GPDs given by the 
Shuvaev transform are not applicable, that is, are in contradiction with the DVCS data. Note, 
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however, that the integral corresponding to the Shuvaev transform is determined directly from 
global diagonal partons, which have a more complicated x behaviour. It is well known that 
already diagonal 'global' data cannot be described by a single Regge pole ansatz. Moreover 
global diagonal partons are poorly described at LO [29]. The NLO corrections are large, since, 
in comparison with LO, they contain qualitatively new 1/z singularities in the splitting and 
coefficient functions. Furthermore, note that the analytic formula for the ratio] (USD, which 
indeed is written for one Regge pole, was given just to qualitatively illustrate the discussion. 

Going to NLO, we see no contradiction of the integral, parameter-free, Shuvaev transform 
approach with the available low x DVCS data [25] . Moreover, the ratios r of skewed-to-diagonal 
PDFs given in Fig. 7 of [M] for the NLO MS scheme are very close to those that we obtain 
using the Shuvaev transform with recent global partons. Table [1] shows the agreement between 
the NLO MS GPDs of |34j and our values obtained using the Shuvaev transform at x — 0.001 
for Q 2 = 10 and 50 GeV 2 . It is clear from the Table that the difference between the GPDs 
based on the Shuvaev transform and those obtained in [31] by the NLO MS fit to the 'skewed' 
DVCS data, is much less than the DVCS error bars. 





Q 2 = 10 GeV 2 


Q2 = 5 o GeV 2 


r G 




r G 


rqs 


CTEQ6.6M 


1.04 


1.67 


1.06 


1.71 


MSTW2008NLO 


1.03 


1.65 


1.05 


1.70 


NLO MS GPDs of [M] 


1.06 


1.76 


1.05 


1.77 



Table 1: Comparison of the skewed-to-diagonal ratio for gluons tq = H g (x,x) / H g (x,0) and 
u, d, s quark singlet rqs = H^(x, x) / xH^(x, 0) at x = 0.001 for scales Q 2 = 10 and 50 GeV 2 . 
Diagonal CTEQ6.6M and MSTW2008NLO partons are used as input to calculate the skewed 
distributions, and the resulting ratios are compared to the ratios as displayed in Fig. 7 of [3J] 
using their NLO MS (£-PW) GPDs. 



6 Conclusions 

We have shown how the skewed (or generalised) parton distributions, Hi(x,^), that are needed 
to describe diffractive processes, can be obtained, for small x, directly from the conventional 
integrated global parton densities, fi(x). This parameter- free method is based on the Shuvaev 
transform, and is applicable to accuracy^ 0(£) at NLO, which should be sufficient for the 
description of all diffractive processes of interest. For the reasons given in the previous section 
the alternative approach of Ref. [31] does not have the advantage of being parameter free. 

5 It is misleading to quote the ratio as "conformal" (as was done for eq.(28) of [34]) since it relies on a single 
Regge pole ansatz, as well as on conformal symmetry. 
6 The accuracy is 0(£ 2 ) at LO. 
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First we identified the reason for the doubts which had been raised concerning the use of the 
Shuvaev transform. Further, we noted that the Shuvaev transform follows from the physically- 
motivated conjecture that the small x input is specified by Regge physics, and so does not 
generate any singularities in the right-half (j > 1) plane in the space-like (x > £) domain. The 
transform therefore has practical applicability. We then investigated the kinematic range of 
reliability of the simplified analytic formula ffTol) . This formula was derived for x — £, assuming 
that the unintegrated parton distributions xf) have a power-like x~ Xi behaviour. The 'analytic' 
formula was only introduced to illustrate some of the qualitative features of GPDs. 

In order to obtain accurate generalised parton distributions at NLO at any small x, £, we 
performed a detailed computation of the integrals of the Shuvaev transform. Finally, these 
computations allow us to provide a readily accessible package which allows the evaluation of 
GPDs for arbitrary small x, £ in the space-like (x > £) domain. These NLO GPDs should 
facilitate the theoretical description of all diffractive processes of interest. 

Appendix: The inverse transform 

Although we do not need the inverse transformation for the calculation of the skewed distri- 
butions for small £ from diagonal partons using the Shuvaev transform, it is still worth noting 
that the inverse problem is solvable. Here we shall briefly clarify the inverse transformation 
and show that the regions \x\ > |£| and \x\ < |£| transform separately. 

Our aim is to find the function f(x) whose Mellin moments, 

M N = [ -x N f(x), 
J x 

are equal to Gegenbauer moments of the GPD H(x,£), 

Gn{0 = 2^+3/2) / d ^4(*/0#M), 

Gn{0 — Mpf. Here for definiteness we consider the non-singlet channel, the singlet one can 
be treated in the same manner. 

With the generating function for Gegenbauer polynomials, 

3 1 f) N 

R{t,x) = 1 - 2tx + t 2 £ 2 , 

the moments take the form 
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where the contour of the t-integration goes in the complex plane around t — 0, while the integral 
over v yields the normalisation. Integrating by parts over y (this is needed to make the kernel 
less singular) we rewrite this expression as 

The contour of t-integration can be shrunk to the cut of the radical in the integrand, 
1 / dt 1 i ,1 , 1 f' 1 dt l-i x 

*2 

where the values ti^ = 2/(^ 2 v)[y ± \Jy 2 — i 2 } determine the real (for \y\ > |£|) interval, in 
which 

i? c (t,y) = uyt - -i 2 v 2 t 2 - 1 > 0. 
These relations show that the desired function is given by the integral transformation 

f(x) = j\yK(x,y)[^- + 5(y - 1) - 5{y + l)]H(y,£), 

with the kernel 

1 /*! a 

K(x,y) = --J ™(l-v)-^[x 2 Y-h(Y)], 

Y = vxy -^ 2 v 2 - x 2 , |y| > iCl- 
Note that for |£| <C 1 the transformation becomes close to the identity, f(x) = H(x) + 0(£). 

If \y\ < |£| the branch points t\^ are complex. In this case the function / will be defined 
on the segment lying in the complex plane. Thus the inverse transform does not exist as a real 
function for the whole interval. 
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